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I. INTRODUCTION 

The thermal conductivity is a property that can provide significant information about 
the nature of superconductors. Even in the early studies of high-T c superconductors 
before good single crystal samples were available, thermal conductivity results indicated 
the anomalous behavior of the normal state transport properties of these oxides and many 
properties of their superconducting state. (See Uher [1] for a review of early thermal 
conductivity studies in the high-Tc superconductors.) Recent low temperature thermal 
conductivity measurements in the normal state of an overdoped cuprate superconductor 
Tl 2 Ba2Cu06+<5 [2] have demonstrated that the fermions which carry heat also carry charge 
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(i.e. the Wiedemann-Franz law is fully obeyed, thus providing no evidence for spin-charge 
separation in overdoped high-T c compounds). In contrast, this same group has observed 
a breakdown of the Wiedemann-Franz law at low temperatures (below T ss 0.15 K) in a 
different high-T c cuprate (Pr 2 _ x Ce :c Cu04_ ( 5) near the optimal doping [3] which suggests 
that superconductivity in the underdoped region of the phase diagram may indeed be 
the result of charged bosons rather than the Cooper pairing most likely responsible for 
superconductivity in the overdoped regime. It was pointed out in Ref. [4], however, that 
the experiments of Hill et al. [3] leave room for possibility that charge- carrying excitations 
are fermionic, but that a subtle transition occurs at T . In any case, such experiments 
help to indicate the value of thermal conductivity measurements to probe the fundamental 
properties of new materials. 

The focus of this work is thermal conductivity measurements in systems that appear 
to exhibit microscopic coexistence of magnetic and superconducting order, two types of 
ordering that are generally antagonistic. Although the traditional transport coefficients, 
Hall effect, thermopower and resistivity, provide little or no information in the super- 
conducting phase below T c , thermal conductivity measurements can still probe both the 
phonons and, at least for a reasonable range of temperature below Tc, the electron quasi- 
particles which are the two primary heat carriers in a superconductor. This may be 
particularly useful for magnetic superconductors in the situation where the magnetic or- 
dering temperature T M is less than the superconducting temperature T c since the electron 
contribution and, perhaps to a lesser extent, the phonon contribution, may both be in- 
fluenced by the magnetic ordering. At high temperature the thermal conductivity also 
provides information regarding the scattering processes for these carriers and, as will be 
illustrated, can provide indications of subtle high temperature phase transitions. Thus, 
thermal conductivity can be an important tool for study of new materials over the entire 
temperature range. 

In the 1980's, discovery of the microscopic coexistence of both ferromagnetic and an- 
tiferromagnetic order in the presence of superconducting order in various Chevrel phase 
and rare earth rhodium boride compounds led to a new view of these two antagonistic 
phases. The experiments generated great excitement, and a large number of theoretical 
predictions of new types of order involving spatial modulation of the two very different, 
mutually exclusive order parameters appeared. Difficulty in synthesis of high quality 
samples of these magnetic superconductors coupled with the excitement of the discovery 
of superconductivity above the boiling point of liquid nitrogen in high-Tc oxide super- 
conductors led to greatly reduced activity in this field. (See Maple [5], Bulaevskii et 
al. [6] and Fischer [7] for excellent reviews of the state of experiment and theory in 
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that era.) The discovery [8-11] of the family of rare earth nickel borocarbide (RNi 2 B 2 C) 
superconductors, which exhibit such a wide variety of interesting phenomena, and the 
widespread availability of good single crystal samples [12] has led to renewed interest 
in questions regarding the coexistence of superconductivity and magnetic order. This 
family exhibits superconducting order only (R = Y, Lu), heavy fermion behavior (R = 
Yb), coexistence of superconductivity and magnetic order (R = Tm, Er, Ho, Dy) and 
magnetic order only (R = Tb, Gd) within the same crystal structure. Although we 
have not found any thermal conductivity measurements in the literature for the Chevrel 
and rare earth rhodium boride magnetic superconductors, perhaps due to the difficulty 
of sample preparation, rather complete data now exists for single crystals of the rare 
earth nickel borocarbides. There are indications of coexistence of magnetic order and 
superconductivity in a new system, the 2122 phase of rare earth rutheno-cuprates, also. 
So far there appears to be only one report [13] of thermal conductivity, in an oxygen 
rich Eui.5Ceo.5RuSr 2 Cu20io-5 compound, but the results of this study suggest strongly 
coupled weak ferromagnet-antiferromagnet/superconductor-normal metal transitions at a 
remarkably high temperature (~ 45 K). 

A brief discussion of the experimental considerations is given in section II with a review 
of the theoretical aspects given in section III. Measurements of the thermal conductivity, 
together with some illustrative data from thermopower and resistivity measurements, for 
the borocarbides and the rutheno-cuprate are discussed in section IV. Section V presents a 
short conclusion. The reader reasonably familiar with thermal conductivity and interested 
primarily in the results for magnetic superconductors may go directly to section IV. 

II. EXPERIMENTAL TECHNIQUES 

We will consider in the following, mainly, the thermal conductivity measurements made 
in the Low Temperature Laboratory at Texas A&M University. For this reason, we will 
describe below in sufficient detail the technique employed in this lab. This technique 
is, however, essentially the same as that used in many other labs for low temperature 
measurements of the thermal conductivity. A more general and thorough account of 
available methods of thermal conductivity measurements can be found in Refs. [14, 15] 

The thermal conductivity of samples was measured by the steady state linear heat flow 
method. One end of the sample was thermally isolated (a known, constant heat input via 
a resistive heater was applied to that end), and the resultant steady state temperature 
gradient was measured. The temperature difference, AT, between the two ends of the 
sample (which was typically less than 5% of the absolute temperature) was measured with 
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a Au-0.007 at. %Fe - Chromel thermocouple. Verification of the temperature difference 
was provided by a Chromel- Const ant an thermocouple for most of the temperature range. 
The temperature differences for each sample and the calibration samples (silver foils) were 
chosen to be approximately the same in the different temperature ranges. The sample was 
completely surrounded by a thermal shield held at the temperature of the cold end. Even 
though heat leaks were minimized by design, the actual heat leaks were measured and 
accounted for in the thermal conductivity data. All this suggests an absolute accuracy 
for the thermal conductivity of ±3 % with the relative precision of ±3 % on a specific 
measurement run. The value of the thermal conductivity k(T) is calculated using the 
relation 

<T) = L±, (!) 

where L and A are length and cross-sectional area of the sample, Q is the power input. 
A detailed discussion of the thermal-conductivity measurement technique employed can 
be found in [16]. 

The above-indicated instrumental error in the thermal conductivity measurements 
should be considered as the optimum or, better to say, minimum one. Some reserva- 
tions, therefore, should be added to this. First, it is usually more difficult to make a 
good thermal contact to an oxide sample than to a metallic one. A somewhat greater 
inaccuracy can be expected in the former case. Second, all the thermal conductance mea- 
surements were made in a vacuum of 1 x 1CT 6 Torr or better, so the convective heat losses 
are negligibly small. The radiative heat losses, due to the small temperature difference 
between the base T-end and the AT-end of the sample mentioned above, should be small 
except perhaps at higher temperatures. As was pointed out in Refs. [1, 14], this radiation 
loss can lead to an overestimate of the measured thermal conductivity above T ps 150 K, 
and the error grows with increasing temperature. Measurement of the actual heat leaks 
without sample, as in the present case, provides partial correction for this error. Lastly, 
it is seen from Eq. (1) that an appreciable error in the absolute determination of the 
thermal conductivity (as well as in supporting measurements of electrical resistivity) is 
due to inaccuracies in the sample dimension measurements. Taking all this into account, 
it can be said that the absolute accuracy for borocarbide and rutheno-cuprate samples 
can be worse than ±3 %, primarily due to inaccuracies in determination of the geometric 
factors of the somewhat irregular single crystal samples and contact resistance for the 
pressed powder rutheno-cuprate samples. 

The electrical resistivity, p(T), was measured for all of the samples studied (except 
the rutheno-cuprate) using the standard four-terminal method. These measurements are 
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helpful and even imperative in many respects. First of all, it gives an indication of 
the quality of the sample. Second, p(T) curves have in general some peculiarities at 
magnetic transitions and show the resistive transition at the superconducting transition 
temperature, T c . Last but not least, p(T) data allow use of the Wiedemann-Franz law to 
estimate the fraction of the total thermal conductivity, that is due to the charge-carrier 
contribution, and, hence, the remainder can thought to be due to phonons and magnetic 
effects, if any. In the case of rare-earth nickel borocarbides, both the resistivity and 
thermal conductivity were measured in the afe-plane. 

III. THEORETICAL CONCEPTS 

For consideration and discussion of available thermal conductivity data for the mag- 
netic superconductors, we will refer to the well established physical mechanisms for heat 
transfer in normal metals and superconductors, outlined in well known books or papers 
[1, 14, 17-27]. With this we will try to understand and explain what is observed in mag- 
netic superconductors and make some conclusions about the properties of heat carriers 
in them. Since the main heat carriers in conducting solids are charge carriers (electrons 
or holes) and lattice thermal excitations (phonons), the thermal conductivity behavior 
reflects properties of these quasiparticles and an interaction between them. In magnetic 
states the magnetic excitations (i.e. magnons) can also participate in the heat transport. 
Hence the thermal conductivity measurements can provide an insight into the nature of 
the superconducting and magnetic states of magnetic superconductors. 

The main (and sometimes rather difficult) task in an analysis of the thermal conduc- 
tivity data for some particular material is to separate the contributions from the two main 
type of the heat carriers - electrons and phonons. This requires, primarily, a knowledge 
about an expected temperature behavior of these contributions (which are quite different 
for electron and phonon heat carriers). Crystal-lattice disorder can strongly effect both 
the electron and phonon parts of thermal conductivity and, therefore, must generally 
be considered. Generally speaking, it is important to find out what scattering mecha- 
nism limits the heat transport for given channel (electrons/phonons) at given conditions, 
which are determined by the phase or state of a sample, its temperature, the magnitude 
of applied magnetic field and other circumstances. 

In the simplest way the thermal conductivity of any type of heat carriers can be 
presented by the equation [14, 17, 19, 20] 

k{T) = ^C v vl = l -C,v\ (2) 
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where C v is the heat capacity per unit volume at constant volume, v is the average 
velocity of the heat carriers, / is the mean free path of the heat carriers (that is the 
average distance they travel between the collisions with any obstacles to heat transport), 
and r is the corresponding relaxation time. Since the quantity D = (l/3)vl is just a 
diffusion coefficient of the heat carriers, it can be said that the thermal conductivity is 
determined by the heat capacity and diffusivity of the heat carriers. In line with this, its 
temperature dependence, k(T), is determined by those of the heat capacity and the mean 
free path of the heat carriers, since v is essentially temperature independent (although it 
can undergo a change at some phase transitions). 

Eq. (2) is applicable for both electron and phonon heat transport in solids. Generally, 
the heat transport via electrons and phonons occurs in parallel. It can be written then 



where /%> and K e present the contributions from the phonons and free charge carriers, 
respectively. Below we consider these contributions separately. The phonon heat trans- 
port takes place in any solid; whereas, the charge carrier contribution depends on the 
carrier density, and is, therefore, negligible in insulators. It should be noted as well that 
different types of relaxation processes, which act simultaneously, can limit the thermal 
and electrical conductivities. As a quite good approximation, Matthiessen's rule [17, 18], 
is applied in this case to describe the combined effect of these processes: 



where r e ff is the effective (or total) relaxation time, and the times r s * correspond to different 
relaxation processes. Since the mean free path is proportional to the relaxation time 
(/ = vt), the same relation can be written for the effective mean free path, l cS . According 
to the Matthiessen's rule, the electrical resistivity can be presented as a sum of the partial 
resistivities, corresponding to different processes of the electron scattering. By analogy, 
it is convenient in many cases to consider the thermal resistivity, W — 1/k, [17, 20] as a 
sum of partial resitivities 



A. Phonon thermal conductivity in non-metallic crystals 

For phonon thermal conductivity, one can use in Eq. (2) the lattice specific heat C p , 
and v may be taken as an average velocity of sound, v s . As a good approximation, the 
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Debye model for the lattice specific heat [17, 19] can be used. According to it, 

C p oc T 3 (T < 6 D ), (6a) 

C p = 3iVA;B (T > 9 D ), (6b) 

where N is the ion density. The main problem is to take into account the most important 
mechanisms of phonon relaxation in crystal solids. In doing so it is convenient to operate 
with the thermal resistivity, which can be presented in Matthiessen's approximation as 

l/«p = W p = W pp + W pi + W pc = -^- 2 (TpjT 1 + T^ 1 + Tp^ 1 ) , (7) 

L/ p f s 

where subscripts pp, pi and pe indicate the phonon-phonon, phonon-imperfection and 
phonon-electron interactions, respectively. 

The phonon-phonon and phonon-imperfection interactions, corresponding to the first 
two terms in Eq. (7), occur in any solid; whereas, the phonon-electron interaction can 
be important only in solids with fairly high charge- carrier density. Let us consider, at 
first, a rather perfect non-metallic crystal, where the main contribution to the thermal 
conductivity comes from the phonon-phonon interaction [14, 17, 19, 20]. These processes 
are those in which two phonons can combine to give a third, and vice versa (so called, 
three-phonon processes). They are restricted by the following selection rules: 

ftuii + T1U2 = hu>3, (8a) 

qi + q2 = q3 + g, (8b) 

where the first condition is conservation of energy at the interaction. In the second 
condition for wave vectors of interacting phonons, g is a reciprocal lattice vector. The case 
g = corresponds to the conservation of momentum (or wave vector). These processes are 
called Normal processes or iV-pro cesses. The case g ^ is determined by the interference 
condition for wave vectors (the lattice acts in that event as a diffraction grating). These 
processes are called Umklapp processes or [/-processes. 

Since iV-processes are energy and momentum conserving they do not contribute to the 
thermal resistivity (see Refs. [14, 17, 19, 20] for thorough explanations). These processes 
can, however, affect the thermal conductivity indirectly, enhancing or modifying the effect 
of other scattering mechanisms, i.e. the phonon-imperfection scattering (see discussion 
of some examples in Refs. [14, 20]). In contrast to iV-processes, the large changes in 
crystal momentum in [/-processes make it possible for them to be an effective source of 
thermal resistivity. The probability of [/-processes is, however, temperature dependent, 
so that they play a dominant part at high temperatures (T > ©d), but are of little 
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importance at low temperatures (T <C ©d), where only iV-processes can occur at an 
appreciable rate. The reason is that at any temperature T, only the so called, thermal 
or dominant phonons with energy Hcu ~ ksT are present in an appreciable number and 
can, therefore, be involved significantly in different interactions with other quasiparticles 
(phonons included). The wave length, A p , and modulus, qT, of the wave vector of thermal 
phonons are given by 

A P = 2^, (9, 

and 

It = -r— = 7T 9°' ( 10 ) 



where go is the Debye wave vector. 

It can be seen from Eqs. (8) that to ensure an [/-process, some of the participating 
phonons should have a large enough wave vector, at least q ~ (l/2)(7 according to [17], 
(i. e., comparable with g D ) and high enough energy (comparable with A; B D ). Thus, the 
[/-processes can have an appreciable rate only at temperatures in the vicinity of or above 

©D- 

Now consider the temperature dependence of the thermal conductivity, k(T), of a 
rather perfect non-metallic crystal. At sufficiently high temperatures (T > © D ), the 
specific heat of the crystal should be constant [Eq. (6b)], and, therefore, the temperature 
dependence k(T) is determined solely by that of the phonon-phonon scattering rate, r^ 1 . 
The dominant contribution to the thermal resistivity in this temperature range is from 
[/-processes. The rate, t~^, should increase with temperature, since the total number of 
phonons (which are scatterers to other phonons which carry heat) is proportional to T 
in this temperature range. It is expected [14, 17, 19, 20] that r^ 1 oc T n with n — 1 for 
the three-phonon processes, which is to say that n p oc l/T. More elaborate theoretical 
calculations [17, 20] give 

Kp oc a ^f*f D , (11) 

1 1 G 

where a 3 gives the volume occupied by one atom, M a is the atomic weight, 7g is Griineisen 
constant. The l/T law can be considered as a good approximation for the lattice thermal 
conductivity of fairly perfect crystals at T > ©d [20]. 

For T < D , the probability of [/-processes drops sharply as temperature decreases. 
In this temperature range theoretical calculations give r pp oc exp(©D/6T) and 

k p ocT x exp(@ D /6T), (12) 

with x and b both of the order of unity [17, 19, 20]. In this case the relaxation time r pp and 
thermal conductivity increase exponentially with decreasing temperature. This increase 
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lasts until the phonon-phonon mean free path, l pp = v s t pp , becomes comparable with that 
of the phonon-imperfection scattering or even with dimensions of the sample (boundary 
scattering). In this case the effective time of phonon relaxation becomes temperature 
independent, and the temperature dependence of k p will be determined by that of the 
specific heat, which reduces with temperature as T 3 for T <C 6d- 

In summary, the k p (T) dependence for fairly perfect non- metallic crystals behaves as 
follows. At high temperatures T > 6 D , it is proportional to 1/T (1/T law). As temper- 
ature decreases below 6d, k p (T) increases exponentially according to Eq. (12) reaching 
a maximum at a temperature, at which the phonon mean free path begins to be deter- 
mined by imperfections or outer boundaries of the sample. With decreasing temperature, 
the mean free path becomes temperature independent, and, therefore, the temperature 
dependence of n p is determined by that of the specific heat, which is proportional to T 3 
for T < D . 

A clear manifestation of this exponential law [Eq. (12)] is confined to range 1/30 < 
T/6 D < 1/10 [20]. This sets an upper temperature limit for the position of the ther- 
mal conductivity peak. The peak can be seen, however, only in rather pure and perfect 
non-metallic crystals. Crystal-lattice imperfections (lattice defects) are effective phonon 
scatterers. They cause an extra thermal resistance which can even suppress the thermal 
conductivity peak completely. Although rather considerable theoretical and experimental 
efforts were made to resolve this problem (see [17, 20] and references therein), the un- 
derstanding level achieved only allows some general discussion and speculation regarding 
analysis of imperfection effects in experimental thermal conductivity data. 

The important feature of lattice defects is whether their linear dimensions are larger 
or smaller than the phonon wavelength. The first type of defects include, for example, 
external or internal (grain) boundaries, dislocations, stacking faults. The second type is 
represented mainly by point defects or small precipitate particles due to phase inhomo- 
geneity. It is obvious that some defects can change their type with changing temperature 
since the wavelength, A p , of dominant phonons depends on temperature as 1/T [Eq. (9)]. 
It can encompass hundreds of interatomic distances for T<6 D ; whereas, it is of the or- 
der of the interatomic distance at T > 6d- This is in sharp contrast with the wavelength 
of electrons in good metals, which is always of the order of the interatomic distance, 
independent of temperature. 

Among other imperfections, the point defects are considered as the most important 
source of thermal resistance. These are impurity atoms, isotopes, lattice vacancies and 
interstitial atoms. The point defects introduce small perturbations of mass, force constant 
and nearest-neighbor distance, causing elastic phonon scattering. Their size is of order 
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of the interatomic distance. The known theoretical results [17, 20] give, in most cases, 
only general ideas about the influence of point defects on the thermal conductivity. At 
low temperature T <C 6d (when A p is well over the size of point defects) the mechanism 
of phonon-point defect relaxation is thought to be similar to that of Rayleigh scattering, 
for which the relaxation rate is proportional to u; 4 [17, 20]. Since uo « k^T jh, for dom- 
inant phonons, the relation, r^ 1 oc T 4 , should hold for the rate of phonon-point defect 
relaxation. In this case, for the low temperature range where phonon specific heat is pro- 
portional to T 3 , the thermal conductivity should be proportional to 1/T. If phonon-point 
defect scattering is dominant at low temperatures, the exponential rise in thermal con- 
ductivity due to [/-processes [Eq. (12)] will be suppressed. At temperatures so low that 
the phonon mean free path, l pi = v s r pi , associated with phonon-point defect scattering, 
becomes comparable with the sample dimensions and, hence, temperature independent, 
the thermal conductivity goes to zero with decreasing temperature according to relation 
k p oc T 3 . 

As temperature increases, the phonon wavelength, A p , decreases, approaching an in- 
teratomic distance. This causes the rate of phonon-point defect relaxation to become 
temperature independent. At high enough temperature this is true for phonon relaxation 
by any lattice defects, regardless of the temperature scattering law for the defects at low 
temperature [20]. It follows from the aforesaid that the phonon-point defect contribution 
to thermal resistivity can be appreciable at high temperatures (T > 9d), where it depends 
only slightly on temperature, while at low enough temperatures the phonon wavelength 
becomes so long that scattering by point defects (and by most of the other lattice defects) 
becomes negligible. It should be recalled that the electron wavelength in good metals 
is always of the order of the interatomic distance, and electron-imperfection scattering 
is temperature independent. For this reason, imperfections remain the only source of 
electron scattering for temperatures so low that the electron-phonon relaxation processes 
are essentially frozen out. This important difference in the behaviors of electrons and 
phonons, as heat carriers, should be taken into account in consideration of the electron 
and phonon contributions to thermal resistance at low temperatures. 

In conductors with a fairly high free electron concentration, the phonon-electron scat- 
tering [see Eq. (7)] can give an appreciable contribution to thermal resistance. This 
problem is, however, so closely connected with electrical and thermal conductivities of 
electrons that it is appropriate to consider it after the main concepts of the electronic 
thermal conductivity will be discussed in the next section. 
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B. Electronic thermal conductivity 



For considerations of electronic thermal conductivity, the general expression given by 
Eq. (2) is used, taking C v to be the electronic specific heat C c and v is the Fermi velocity 
v F , i. e. 

K e (T) = ^C c v F l = ^C c v 2 F r. (13) 
The electronic specific heat in the free-electron model [17-19] is given by 

C e = Y (j^p) n e k B = 7 T, (14) 

where S F is the Fermi energy, n c is electronic density, and 7 is the Sommerfeld coefficient. 
The linear dependence on temperature is quite general, however, and holds for any Fermi 
liquid, independent of the free-electron model. 

It is of fundamental importance that the same scatterers (phonons and lattice defects) 
determine the resistance to both the charge and the heat transport of electrons. The 
electronic conductivity, a, in the frame of free-electron model is given by [17-19] 

n c e 2 r 

a= . 15 

m 

When expressions for electronic thermal and electrical conductivities [Eqs. (13) and 
Eq. (15)] are compared to one another, it is apparent (using Eq. (14) for electronic specific 
heat as well) that K c /(aT) should be constant, if the electron relaxation times (the total 
relaxation time, r) are equal for both kinds of electron transport. This simple relation is 
the famous Wiedemann-Franz (WF) law. In the free-electron model this law is expressed 
as [17-19] 

£ = £0, (16) 

where L = (tt 2 /3) (/cb/c) 2 is called the Lorenz number. 

The WF law forms the basis for separation of the electron and phonon contributions to 
thermal resistance during analysis of experimental data on thermal conductivity. This law 
appears, however, to be true only at the lowest temperatures, where electron-imperfection 
scattering is dominant, and at high temperatures (T/Od ~ 1). In the intermediate 
temperature range, a significant violation of the WF law takes place. The reason is that 
the electrical and thermal conductivities are determined by different types of electron- 
phonon collisions in this intermediate temperature range. To make this point more clear, 
let us write in the Matthiessen's approximation the more detailed expressions for the 
thermal and electrical resistivities: 

l/«e = W C = W cp + W ci = -^ [r^ 1 + t cp(k) -\T)} , (17) 
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P = Pc P + Pci = 2 [V* + M'fW] , (18) 

where subscript ei indicates the electron-imperfection interaction (which is temperature 
independent), and subcripts ep(/t) and ep(cr) indicate electron-phonon interactions, cru- 
cial for thermal and electrical conductivities, respectively. Of course, some other sources 
of electron scattering can be indicated as well (for example, electron-electron scattering 
or the spin disorder in paramagnetic, ferromagnetic or antiferromagnetic states), but the 
two mentioned above are considered as the most significant. 

The times T ep ( CT ) and t cp(k ) are characterized by the following temperature behaviors 
[17-19]: 

[t 5 ifT<e D , 
% CT )( T )°M » (19) 

(T ifT>6 D ; 
IV 3 ifT<6 D , 

{> [t ifT>e D . 

It is seen that temperature dependences of these times are the same near or above 0d, 
but are quite different below it. Consider briefly the sources of this diversity. When 
an electron with wave vector k is scattered by a phonon into the state k', the electron 
either absorbs or emits a phonon of wave vector q. These collisions are determined by 
the following rules [17-19]: 

k' = k±q + g, (21a) 
£ k > = £ k ± hoo^. (21b) 

The first rule maintains that momentum is conserved, up to the addition of any arbirary 
vector of the reciprocal lattice. In the same way as for the phonon-phonon scattering (see 
Sec. Ill A), the cases g = and g ^ correspond to N- and [/-processes, respectively. It 
should be noted that the [/-processes of the electron-phonon interaction are very improb- 
able at low temperature, but they can increase somewhat the electronic electrical and 
thermal resistivities at fairly high temperatures (above ~ 0.2 ©d)- The second rule is the 
requirement for conservation of energy in the interaction. 

The maximum phonon energy is about /cb©d, which is much less than the electron 
energy, £p. For this reason, the electron-phonon collisions could be considered as quasi- 
elastic ones. For electronic thermal conduction, however, it is more important to compare 
the change in electron energy in a collision (A£ q = |£ k — £ k '|) with k^T (since electrons 
in metals can change their energy only in narrow band near the Fermi level with a width 
about kv,T). Following this criterion, the collisions with A£ q ~ k^T are considered 
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as inelastic ones while those with A£ q < /c B T, as elastic [17-19]. These collisions are 
quite effective for electronic thermal resistance. At low temperatures, T « D , they 
are inelastic (since the energy of the dominant phonon is about k^T in this temperature 
range) with a collision rate (or an inverse relaxation time, l/r cp ) proportional to T 3 
[17, 19]. This corresponds to the low temperature behavior of time r ep ( K ) in Eq. (20), 
governing the electronic thermal resistance. 

At low temperatures, T<6 D ,a single collision, which must be an iV-process according 
to the selection of Eq. (21), is, however, ineffective in providing electrical resistance. The 
reason is that the electron direction of motion can be changed only slightly in this N- 
process. The angle between vectors k and k' (the angle of scattering) is 9 ~ (/tAf = 
T/Od, that is, very small. An electron must endure many of such collisions to gain a large 
scattering angle, and the probability of that event is proportional to (T/6d) 2 [17-19]. For 
this reason, the effective rate of the electron-phonon relaxation for electrical resistance is 
t~J ct) (T) oc T 3 (T/0 D ) 2 oc T 5 , which was already indicated in Eq. (19). 

At high temperatures, T > D , any electron-phonon collision is large- angle and elastic. 
Since the total number of phonons is proportional to T at this temperature range, the 
relation oc T holds for both electron charge and heat transport. 

We can survey now a temperature dependence of the electronic thermal conductivity, 
taking into account Eqs. (13), (14) and (17). At temperatures so low that the elastic 
and temperature independent electron- imperfection scattering is dominant, the thermal 
conductivity behaves as k c (T) oc T. In this temperature range the WF law [Eq. (16)] 
holds. For higher temperatures, the electron-phonon scattering becomes dominant. At 
T <C Od, this scattering is inelastic with rate, T e ~^ oc T 3 , relevant for heat transport, and 
the rate t~}^ oc T 5 , relevant for charge transport. This leads to the relation k c (T) oc T~ 2 . 
Due to different electron-phonon relaxation rates, relevant for charge and heat transport, 
the WF law is not obeyed in this temperature range, so that the ratio k c /(<jT) is not 
constant, but is proportional to T 2 . At high temperatures, T > D , electron-phonon 
collisions are elastic with rate oc T, which gives relations, p oc T, for electrical 
resistivity and, k c (T) = const, for thermal conductivity. The WF law holds in this 
temperature range. For fairly perfect metals with high electron density, the electron 
contribution dominates in thermal conductivity, so that the outlined features of n e {T) 
dependence are immediately evident from measured n(T) curves. An example of such 
behavior is presented in Fig. 1 for a well annealed silver foil [16]. 

It is seen that a concurrence of the electron-imperfection and electron-phonon interac- 
tions gives rise to a maximum (a peak) in the K e {T) dependence. This peak can be very 
sharp for fairly perfect and pure metals, where the peak position, T max , is usually below 
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FIG. 1: k c (T) dependence of annealed silver foil. 

0.1 ©d- Lattice defects induce the following: (i) an increase in slope of the linear part of 
K e (T) dependence; (ii) an increase in T max ; and (iii) a decrease in the peak height up to a 
total smearing of it for strong enough lattice disorder. Some clear demonstrations of this 
imperfection effect for typical metals can be found in [14, 15, 18, 20, 28]. It is apparent 
also [see Eq. (17)] that in "bad" metals with strong electron- imperfection scattering the 
WF law is obeyed better (even in the intermediate temperature range where inelastic 
electron-phonon scatering occurs) than in pure metals with a perfect crystal lattice. 

C. Effects of phonon-electron interactions in thermal conductivity 

The phonon-electron scattering can limit the phonon thermal conductivity (this point 
was briefly mentioned at the end of Sec. Ill A). This interaction leads to absorption 
or emission of phonons by electrons and is restricted by the above-mentioned rules for 
electron-phonon interaction [Eq. (21)]. The thermal resistivity, determined by these pro- 
cesses, is labelled W pe in the general expression (7) for phonon thermal resistivity with 
an associated relaxation time r pe . The known theoretical examinations [14, 17, 20] have 
shown that tZ}(T) oc T if T <C 6d, and is constant if T > D . This leads [taking into 
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account Eqs. (6) and (7)] to 




T -2 



const if T > D . 



if T<6 D , 



(22) 



Equations (22) actually results from the more general relation [17] 
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where n a is the number of conduction electrons per atom. Besides the electron-phonon 
iV-processes, this equation takes into account the [/-processes as well, and that causes an 
appearance of the nonequality sign. It is easy to verify that in the case of the equality 
sign (that is, ignoring the ^/-processes) the relationships (22) can be derived from Eq. (23) 
for low (T <C D ) as well as for high (T > D ) temperatures. In the former case, the 
relationships p oc T 5 and C p oc T 3 should be taken, and in the latter case, p oc T and 
C p = 3iVA;B. 

Equation (23) has a doubtless similarity with the WF law presented by Eq. (16). At 
high temperatures (C p = 3NIcb) it is possible to rewrite it as 



The right-hand side of the equation is quite close to the Lorenz number, L = 
(n 2 /3)(kB/e) 2 . This implies that the phonon thermal conduction at high temperatures is 
comparable with the electronic thermal conduction, if both of them are determined by 
the electron-phonon scattering only [17]. 

D. Superconductivity effects in thermal conductivity 

In superconductors, as temperature crosses T c from above, some number of electrons, 
n s , becomes superconducting, while the rest of them, in number n n , would remain as before 
in the normal state. The fraction of superconducting electrons [given by x = n s / (n s + n n )} 
increases continuously with decreasing temperature from x = at T = T c to x = 1 at 
T = 0. The superconducting electrons are in a bound state (Cooper pairs), in which they 
are unable to transport entropy or interact with phonons. But the rest of the electrons 
(the fraction, 1—x) being normal, remain heat carriers and can interact with the phonons. 

Even in the frame of this simple picture (which is in the spirit of the two-fluid model of 
Gorter and Casimir) some quite definite predictions about the effects of superconductivity 




(24) 
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on the thermal conductivity can be made. Before doing this, let us denote the electronic 
thermal conductivity in the superconducting state as «|, and that in the normal state 
(induced, for example, by high enough magnetic field) as The same kind of notations 
(ftp and ftp 1 ) will be used in the following for the phonon part of the thermal conductivity 
or even for the total thermal conductivity (k s and K n ). It is clear that the ratio k|//c£ 
must reduce continuously with decreasing temperature below Tc from the initial value 
(equal to unity) at T = Tq to values much less than unity at low enough temperature. All 
this takes place due to a decrease in the fraction of the normal electrons when going below 
Tq. The same reason (that is a decrease in the number of the normal electrons below T c ) 
can, however, induce an increase in the ratio, ftp/ftp 1 , for the phonon contribution, causing 
it grow far above unity with decreasing temperature in the range not too far below T c . 
In this case a reduction in the number of normal electrons (which are the main phonon 
scatterers in metals for low temperatures) can cause a decrease in the phonon-electron 
relaxation rate, l/r pe , and, hence, an increase in the phonon thermal conductivity (it 
follows from Eq. (23) that K pc oc l/n e ). 

The basic theoretical results on the superconductivity effects in the thermal conductiv- 
ity were obtained mainly in the frame of the BCS model [21-27]. Consider them separately 
for the electron and phonon contributions to the thermal conductivity. It turns out to 
be rather important for the electron heat transport whether the critical temperature Tc 
falls below or above the peak in the temperature dependence of K e (T). In the former 
case, the normal electrons (or as it is often said, quasiparticle excitations) are scattered 
predominantly by the lattice imperfections, while in the latter case, by phonons. For both 
cases rather cumbersome expressions for temperature dependences of k b J k£ were derived 
which can be found in Refs. [1, 21-27]. We shall restrict ourselves to consideration of 
the graphical representations of these dependences shown in Fig. 2(a). When defects are 
the dominant electron scatterers, the function k%/ ' = f(T) has a zero slope at T c [solid 
curve in Fig. 2(a)]; whereas, in the case that electrons are scattered mainly by phonons 
[dashed curve in Fig. 2(a)], the slope at T c is rather large (about 1.62 according to Ref. 
[24]). Theories [21-27], including the so called BRT model [22, 24, 26], agree well with 
experiment for the case of predominantly electron-imperfection scattering that is found 
for pure low-T c superconductors with high electron density like Al, In, Sn and some others 
(see [14, 24, 27]). For the case where elecron-phonon scattering dominates, the agreement 
between theory and experiment is not so convincing. 

It should be noted that for a comparison between the theories [21-25] and experimental 
data one must know exactly the "normal" behavior of thermal conductivity below T c 
since a theoretical expression for — f(T) is used for the comparison. For a weak 
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superconductor with low Tc, the "normal" thermal conductivity can be measured below 
Tq under a magnetic field which is high enough to suppress superconductivity. This is 
not feasible, however, for superconductors (like cuprates) with high Tq and huge critical 
magnetic fields. Besides, when T c is high enough (as in the cuprates) and falls into the 
temperature range where electron-phonon scattering dominates, the WF law fails and 
does not give a good estimate of the electronic part of the thermal conductivity from the 
total measured thermal conductivity (see Sec. Ill B). The comparison with BRT and 
other models in the case of superconductors with high enough T c poses, therefore, great 
difficulties. 




FIG. 2: Left panel (a) presents theoretical temperature dependences of the ratio of electronic 
thermal conductivities in superconducting and normal states, k%/k^, below Tq for the cases 
when the electrons are scattered predominantly by defects (solid line) or by phonons (dashed 
line) [21-25]. Right panel (b) presents schematic sketches of temperature dependences of the 
total thermal conductivity below Tq for different relationships between the electron and phonon 
contributions to the thermal conductivity. Expected behaviors for the cases when the electrons 
or phonons dominate the thermal conductivity are presented by curves 1 (doted line) and 3 
(dashed line), respectively. Curve 2 (solid line) presents the intermediate case when both kinds 
of the heat carriers compete on equal terms in the thermal conductivity. 

The phonon contribution to heat transport can be quite appreciable for some materials 
and be comparable or even higher than that of the electrons. In this case the thermal con- 
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ductivity behavior below Tq can be partially or totally determined by phonon transport. 
Figure 2(b) presents schematic sketches of possible temperature dependences of the total 
thermal conductivity, /t s , below Tq for different relationships between the electron and 
phonon contributions to the thermal conductivity. A behavior of k s {T) for a limiting case 
when the electron contribution dominates is presented by curve 1 (doted line). This type 
of behavior is expected for rather perfect superconducting metals with a high electron 
concentration. Another limiting case, when phonons dominate in the thermal conductiv- 
ity, is illustrated by curve 3 (dashed line). In this case (as was already mentioned above), 
the phonon thermal conductivity increases below T c with decreasing temperature due 
to the decrease in the number of normal electrons with a corresponding decrease in the 
phonon-electron relaxation rate. At low enough temperature the phonons begin to be 
scattered primarily by crystal boundaries, so that the phonon mean free path becomes 
temperature independent. After that, according to general relation (2), k s (T) decreases 
to zero with decreasing temperature with the resulting formation of a maximum in n s (T) 
dependence [curve 3 in Fig. 2(b)]. This type of the K s (T)-behavior is expected in dis- 
ordered metals and alloys with high electron-imperfection scattering (that depresses the 
electron contribution to the thermal conductivity), and was seen in some disordered low- 
T c superconductors [25, 29]. This behavior was found to be in good agreement with the 
theory of Geilikman and Kresin [25]. The phonon thermal conductivity may be dominant 
also in rather perfect metals which have a low free electron density. This is characteris- 
tic of high-Tc cuprates where K s (T)-behavior, similar to that of curve 3 in Fig. 2(b), is 
frequently observed [1]. 

The curves 1 and 3 in Fig. 2(b) present the limiting cases when one type of heat 
carriers (electrons or phonons) dominates the thermal conductivity. It is apparent that, 
for many materials, the electron and phonon contributions will be of the same order. This 
will lead to dependences which are intermediate between the two limiting cases. Curve 2 
(solid line) in Fig. 2(b) presents one of the possible intermediate cases when both kinds 
of the heat carriers compete on equal terms in the thermal conductivity. Some examples 
of this kind of k s (T) behavior for low-T c superconductors can be found in Ref. [25]. 

IV. RESULTS AND DISCUSSION FOR MAGNETIC SUPERCONDUCTORS 
A. Rare-earth Nickel Borocarbides 

In this section we consider the thermal conductivity data for rare-earth nickel borocar- 
bides, RNi 2 B 2 C, with R = Y, Lu, Yb, Tm, Er, Ho, Dy, Tb or Gd (which will be referred 
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below to as borocarbides). Although only those members with R = Tm, Er, Ho and Dy 
are magnetic superconductors, for the sake of completeness and to illustrate similarities 
for the members that exhibit only one of these behaviors (only superconductivity, only 
magnetism) or even neither, all members will be discussed. Before going directly to this 
matter, it is appropriate first to briefly review the basic properties of these compounds. 
For extended reviews of experimental and theoretical findings for borocarbides see Refs. 
[30-33]. 

1. General Properties 

The crystal structure of RNi 2 B 2 C is a body-centered tetragonal with space group 
IA/mmm [34]. It is a layered structure in which Ni 2 B 2 layers are separated by R-C 
planes stacked along the c-axis. The Ni atoms are tetrahedrally coordinated to four B 
atoms. Despite the layered structure, the electronic properties of borocarbides appear to 
be essentially three-dimensional according to electronic band structure calculations [35], 
supported by the spectroscopic studies [36] . These compounds are rather good metals with 
a large electron density of states N(£ F ) at the Fermi level. The dominant contribution 
to N(£ F ) comes from Ni 3d-states, but some contributions come from all other atoms as 
well [35]. Resistivity and magnetoresistance of RNi 2 B 2 C borocarbides (R = Y, Lu, Er, 
Ho) [37] have not revealed appreciable anisotropy between measurements normal to the 
afe-plane (along the c-axis) and in-plane, consistent with the three-dimensional character 
of the electronic properties. Small deviations from isotropic behavior were observed below 
150 K for the compounds containing magnetic ions (Er and Ho), which was attributed in 
Ref. [37] to crystal-electric field effects. 

Among the RNi 2 B 2 C family, the compounds with R = Y, Lu, Yb are non-magnetic 
(the last is a heavy-fermion system), but all others (R = Gd, Tb, Dy, Ho, Er, and 
Tm) are magnetic. Their magnetic properties are determined by the localized electrons 
in the 4/-shell of R atoms. Long range magnetic order is due to the indirect RKKY 
interaction, mediated through the conduction electrons. This gives rise to different types 
of antiferromagnetic ( AFM) order of the 4/-ions [38] . Borocarbides with R = Tm, Er, Ho, 
Dy show coexistence of superconductivity and magnetic order. More information about 
superconducting and magnetic properties of these compounds can be found in Table I. 

A Cooper pair consists of two electrons with equal and opposite moment, and with 
opposite spins (the total momentum and spin of the pair are zero). Any perturbation 
which acts with opposite signs (or with opposing force) on the two members of a Cooper 
pair can destroy this pair, producing a, so called, pair breaking effect. For example, 
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TABLE I: Low temperature phase states of the RN12B2C bulk compounds (taken from Refs. 
[31, 33, 38] unless otherwise referenced) . R = Y or rare-earth element, Tq is the superconducting 
transition temperature, Tm is the magnetic transition temperature, Tk is Kondo temperature. 
Column "Magnetic Transition" indicates magnetic states into which the compounds can be 
transformed (AFM = antiferromagnet, MAFM = modulated incommensurate antiferromagnet 
and WFM = weak ferromagnet) , q is the wave vector of the modulated magnetization (in 
units of reciprocal lattice parameters) and the symbol ||| indicates magnetic moment directions. 
Abbreviation HFS denotes heavy-fermion system. 
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an external or the internal (in magnetic materials) magnetic field exerts an orbital pair- 
breaking effect since the field acts with opposing force on the two electron momenta in 
the pair. If some ions in a metal system have a magnetic moment, ion spins will act with 
opposite sign on the electron spins in the pair (the magnetic pair-breaking). Generally 
speaking [42] , such types of perturbations lead to breaking of the time-reversal properties 
of the system and, hence, can cause a strong decrease in T c or even total suppression of 
superconductivity. 

In the case of rare-earth compounds, localized 4/-electrons should undoubtely exert 
a pair-breaking effect on the superconductivity. Since, however, the 4/-electrons are 
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strongly localized in deep inner 4/-orbitals, their interaction with conducting electrons 
can be rather weak, thus, permitting coexistence of superconductivity and the long range 
magnetic order. Such a situation is believed to take place for rare-earth rhodium boride 
compounds, Chevrel phases [5-7], and borocarbides [30]. Coexistence of superconductiv- 
ity and AFM order was justified theoretically rather long ago [43]. The coexistence is 
possible if the AFM exchange field averages to zero within the superconducting coherence 
length, £. AFM order is, however, expected to have a profound effect on superconducting 
and transport properties, and this was found in the above-mentioned magnetic supercon- 
ductors. 

The theoretical and experimental studies generally indicate that the borocarbides are 
conventional electron-phonon mediated superconductors with s-wave symmetry for the 
order parameter [30, 36]. On the other hand, some members of the borocarbide family 
manifest properties which are suggestive of unconventional or exotic superconductors 
[32, 44-46]. In spite of this, the borocarbides are reasonably considered to be conventional 
superconductors. 

Rare-earth borocarbides are rather good conductors, comparable to transition metals 
or their alloys. They have high electron density. According to Ref. [47], the value of n c 
for Y-, Ho-, and Gd-based borocarbides are 2.63, 3.12 and 5 x 10 22 cm" 3 , respectively. 
These values are only moderately less than those of simple superconducting metals, such 
as Al, Pb or Sn [19]. The Fermi velocity values, 3.6 x 10 7 cm/s (calculated for LuNi 2 B 2 C 
by Pickett and Singh [35]), and 4.2 x 10 7 cm/s (found by de Haas- van Alphen studies 
in YNi 2 B 2 C [48]) are, however, clearly much less than the typical value vp ~ 2 x 10 s 
cm/s for simple metals [19]. Measurements of transport properties of borocarbide single 
crystals [31, 37, 49] have shown that the resistivity at room temperature, prt, is typically 
in the range 35-70 pflcm, and that the residual resistivity, po, is in the range 2-5 pfl cm 
(T < 10 — 15 K), so that the ratio par/po is in the range 10-30. These conducting 
properties, though rather good, are not so good as those of the pure simple metals, for 
which po can be as low as 0.01 pQcm and Prt/po can be as high as a thousand [50]. 
It thus follows that the relative electron contribution to the total thermal conductivity 
in borocarbides should be much less than that in pure simple metals. For this reason, 
some manifestations of the phonon contribution in the behavior of the total thermal 
conductivity of borocarbides is expected. In particurlar, due to the rather high residual 
resistivity, p , the phonon effects may show themselves below T c . The high p - values 
even in single-crystals are possibly determined by some types of inhomogeneities and 
point defects. The known studies [51] suggest vacancies at the boron/carbon sites, which 
have no appreciable influence on bulk superconducting properties. The vacancies can raise, 
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however, p -values considerably, depressing in this way the electronic thermal conductivity 
and enhancing the role of the phonon thermal conductivity in the low temperature range. 

The temperature dependence of the resistivity in borocarbides deserves some attention 
since it may be important for consideration of their thermal conductivity The experi- 
mental p(T) curves are found to be approximately linear in the range 100-300 K (see p(T) 
plot for LuNi 2 B 2 C in Fig. 17) [31, 49]. Below 50 K, p(T) is essentially non-linear and, in 
a rather narrow range 1.25 T c < T < (30-40) K, can be approximated by p(T) = p + aT p , 
where p- values for RNi 2 B 2 C compounds are in the range 2.0 < p < 2.6 [31, 49]. 

Since the specific heat enters the general expression (2) for the thermal conductivity, 
a brief mention should be made of this point. Most studies focussed on the specific heat 
in the neighborhood of T c (above and below it) [30, 52-54]. The lattice specific heat, 
C p , in fairly good metals begins to exceed the electronic contribution at a temperature 
which is a few percent of the Debye temperature [19], and the same must be true for the 
borocarbides as well. In the Debye model, C p tends to a constant value when T approaches 
6d- It is essential, therefore, to know the Go-values for the compounds considered. In 
the specific heat measurements the following values of 6 D were obtained: 345 K for R 
= Lu [52], « 540 K [53] or « 490 K [54] for R = Y, and 320 K for R=Tm [54]. It is 
evident that On-value for other borocarbides fall within the outlined temperature range. 
Since most of the thermal conductivity measurements considered below were made in the 
range 4.2-300 K, it is important to know how C p behaves in this range, especially in the 
high-temperature range 100-300 K. Judging from the above-indicated On-values, a fairly 
large rise in C p with increasing temperature should take place in the high-temperature 
range. This was found in LuNi 2 B 2 C [55], where C p has nearly doubled in the range 
150-300 K. Although no other specific-heat studies were done (to our knowledge) in the 
high-temperature range, the same kind of the C P (T) behavior must be expected for the 
other borocarbides as well. 

2. Thermal Conductivity 

The thermal conductivity studies of RNi 2 B 2 C compounds are few in number [56-62]. 
Some of them [56-58, 60, 62] have been made only at very low temperature. Below we 
consider results on single-crystal borocarbides with R = Y, Lu, Yb, Tm, Er, Ho, Dy, Tb, 
and Gd. Certain of the results have been partially presented in Refs. [16, 59, 61]. The 
samples were provided by P. C. Canfield and his research group at Iowa State University. 
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a) Borocarbide magnetic superconductors 

Let us start with the magnetic superconductors (R= Dy, Ho, Er, Tm). The tempera- 
ture behavior of their thermal conductivity is presented in Figs. 3-10. Half of the figures 
present expanded views of k(T) and p(T) (and even the absolute thermopower, S(T), in 
Fig. 4) dependences at low temperature. This allows one to follow closely the changes 
in these properties at phase transitions. On the panels with p(T) curves, the values of 
Prt/Po are indicated for characterization of the sample quality. 3 On all figures, /t to tai 
stands for the measured thermal conductivity, n e is the WF estimate for the electronic 
contribution (based on the measured p for these samples), and the remainder, K tota \ — tz c , 
can be thought as being the phonon contribution. 

The DyNi 2 B 2 C compound transforms from the paramagnetic (PM) into AFM state 
at the Neel temperature (T N m 10.3 K), which is considerably higher than T c ps 6.2 K 
(see Table I and Refs. [38, 63]), so that this compound goes to the superconducting 
state already being in the AFM state. The PM-AFM transition is of first order and 
characterized by very distinct and intensive peak in the temperature dependence of the 
specific heat [63] . It can be seen from Fig. 4 that the transition results in a reduction by 
nearly half in the resistivity, and in an appreciable increase in the thermal conductivity 
and thermopower (in the latter case, the increase is just enormous). The main reason for 
these dramatic changes is a decrease in the rate of electron scattering by the spin disorder 
in response to the PM-AFM transition, thus producing the pronounced structure in k(T), 
p(T) and S{T) shown in Fig. 4. It is known that spin disorder can give a considerable 
contribution to the resistivity above T N (see the discussion in Ref. [31]). Actually, a 
relevant spin-disorder term, p sd (with the electron relaxation rate, r^ 1 ) must be added 
into Matthiessen's relation (18) for T > T N . This term can be reasonable above T N , 
but, for ideal spin alignment below T N , it may be thought to be equal to zero. The 
same reasoning is applicable to the thermal conductivity with the resulting addition of 
the term r^ 1 to the relation (17). The relaxation rate, r^ 1 , is generally considered to be 
temperature independent above Tn [64]. 

It is evident from Fig. 4 that the main change in the K, tota \(T) at the PM-AFM 
transition can be attributed to that of the electronic contribution, K e . Nevertheless, a 
small dip can be distinguished in the temperature curve of the phonon contribution, 



3 po here is the value at the temperature Tq at the onset of the resistive superconducting transition, not 
necessarily the residual resistivity. For R = Tm, Er, Ho it includes the contribution from spin disorder 
scattering. 
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FIG. 6: k(T) (a) and p(T) (b) for HoNi 2 B 2 C at low temperature. 
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FIG. 8: k(T) (a) and p(T) (b) for ErNi 2 B 2 C at low temperature. 
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FIG. 10: k(T) (a) and p(T) (b) for TmNi 2 B 2 C at low temperature. 
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Ktota\—K>e, in the vicinity of Tn, as well. The phonon thermal conductivity can be respon- 
sive to a magnetic transition when this is accompanied by strong enough magnetoelastic 
(magnetostriction) effects. 

There is no easily discernable change in K total (T) at T c (Fig. 4). It is possible that 
rather small changes in slope of k(T) at T c is somewhat hidden on the background of the 
rapid variation in k(T) due to the PM-AFM transition. This variation can be appreciable 
even well below T N down to T c . This is supported by the point that fairly large changes 
of p(T) and S(T) due to this transition occur in the temperature range between T N and 
T c [Fig. 4(b)]. No sharp change should be expected either if the electron contribution, re e , 
dominates, and the electrons are scattered predominantly by defects. No abrupt change 
in slope of K e (T) at T c occurs in that case [see Fig. 2(a) and discussion in Sec. HID]. 

The HoNi 2 B 2 C compound has three magnetic transitions below T c ~ 8.5 K (see 
Table I and Refs. [38, 65, 66]), so that a coexistence of bulk superconductivity and 
local long-range magnetic order takes place. The p(T) shows a sharp superconducting 
transition [Fig. 6(b)], but K total (T) does not show a noticeable change in behavior at T c 
[Fig. 6(a)] (the same is found in Ref. [56]). This may be related to defects as the dominant 
electron scattering mechanism as discussed for the case of DyNi 2 B 2 C (see above), but it 
is also may be connected with an anomalous superconducting state in the temperature 
range between Tc ~ 8.5 K and 6.5 K, which was indicated by point-contact tunneling 
measurements [67]. Of these magnetic transitions, only the transition to the simple AFM 
state at T N 5.2 K produces an appreciable feature in the n tota \(T) curve [Fig. 6(a)]. 
The other two [their transition temperatures are denoted by T M1 and T M2 in Fig. 6(a)] 
have a minimal effect on the thermal conductivity. A change in slope at Tm 2 is clearly 
visible, but any feature at T M1 is within the noise. It is interesting that the specific heat 
(C) behaves somewhat in a similar manner at these transitions [65, 66]. Whereas, a sharp 
and intensive peak in C(T) is found at T = T N [65, 66], there are only weak discernible 
shoulders in the C(T) dependences at the temperatures T M1 and T M2 [30, 65]. The reason 
is that AFM transition at T N is of first order [66]; whereas, the other two are of second one. 
It is clear that first-order magnetic transitions, characterized by discontinuous changes in 
the magnetic entropy (see Refs. [63, 65] for borocarbides) and in the compound density, 
should lead to more significant change in the thermal conductivity of the borocarbides than 
the second-order ones (characterized by discontinuous changes in the second derivatives 
of the Gibbs free energy, such as the volume thermal expansivity or the specific heat). 

The ErNi 2 B 2 C borocarbide transforms with decreasing temperature first into the 
superconducting state at T c ps 11 K, then into a MAFM state at T N xs 6.8 K (first order 
transition), and finally into a WFM state at T WF ss 2.3 K (second order transition) (see 
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Table I and Refs. [38, 40, 68]). For this compound a rather distinct change in the slope of 
Ktotai(T) at T c can be seen [Fig. 8(a)], consistent with measurements on polycrystalline 
ErNi 2 B 2 C sample by Cao et al. [57]. But no features can be distinquished in the K, tot3b \(T) 
curve at T N and T WF . 

The TmNi 2 B 2 C compound is superconducting below Tq ~ 11 K, and undergoes a 
MAFM transition (of first order) below T N « 1.5 K (see Table I and Refs. [38, 54, 69]). A 
rather clear kink is seen in the K totSL i(T) curve at T c , which corresponds in Fig. 10(a) to 
the temperature of the superconducting transition from the p(T) data. A sharp increase 
in the thermal conductivity is readily seen at about 1.4 K, that is slightly below T N [see 
insert in Fig. 10(a)]. It is particurlarly remarkable that fttotai doubles in a very short 
temperature range. It is clear that this increase is determined by the transition to the 
AFM state. T N is, however, much less than T c (T N /T C ~ 0.136) and, therefore, the 
fraction of the normal electrons (which remain heat carriers below T c , see Sec. HID) 
must be very small below 7n- For this reason, the explanation given above in the case 
of DyNi 2 B 2 C (an abrupt decrease in the rate of electron scattering by the spin disorder 
below T N ), is inapplicable for TmNi 2 B 2 C. Really, the theoretical calculations [see Refs. 
[22, 24] and Fig. 2(a)] show that K e is negligibly small for T N /T C < 0.2. 

The transition of TmNi 2 B 2 C into the AFM state does not destroy its global super- 
conductivity in the sense that the resistivity remains zero [54]. But it could suppress the 
superconducting state in such way (causing the effective T c to become less than 11 K) 
that the fraction of the normal electrons would increase rather abruptly below T N . To- 
gether with disappearing of the spin-disorder electron scattering below T N , this could lead 
to an appreciable increase in K e . A considerable decrease in the upper critical field H c2 
(parallel to the ab plane) found in TmNi 2 B 2 C below 3 K [31, 70] gives some support to this 
suggestion. According to Ref. [69], the Tm moments order below 1.5 K with an incom- 
mensurate magnetic structure consisting of Tm moments, aligned ferromagnetically along 
the c axis in the (110) planes, and the magnitude of the moments modulated sinusoidally 
along the diagonal of the ab plane. This magnetic structure should, however, support 
magnons. The formation of the magnons, which are heat carriers and whose number is 
especially large in the temperature range just slightly below Tn, also could explain the 
sharp rise in K to tai below T N . It is hard to judge at the moment, however, the possible 
relative contributions, if any, of these two inferred mechanisms without some additional 
experimental studies. 

Let us survey now the general features of the thermal conductivity behavior in the 
borocarbide magnetic superconductors (R = Dy, Ho, Er, Tm) in the temperature range 
above both, T c and T M (Figs. 3, 5, 7, and 9). It can said that the electron contribution to 
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the thermal conductivity, K e , clearly dominates over the phonon one, K tota i-K c , only at low 
temperature (< 50 K), but with increasing temperature the relative phonon contribution 
grows significantly, so that at room temperature these contributions are quite comparable 
(if not nearly equal). The K e (T) is approximately constant for T > 100 K, as expected 
according to Eq. (13) for electronic thermal conductivity and the WF law (16) in the case 
of linear p(T) dependence (which is approximately obeyed in this temperature range for 
all of the borocarbides mentioned). But there is no peak in the n e (T) dependence which 
is expected below 0.1 6d for good enough metals (see Sec. Ill B). The reason is that the 
borocarbides are maybe not good enough metals to the extent that they should show this 
peak. Really, the highest thermal conductivity magnitude of borocarbides is generally 
below 0.3 W/cmK (see Figs. 3, 5, 7, 9 above and figures for other samples below). At 
the same time, the silver foil presented in Fig. 1 has a peak of about 14 W/cmK at 
T miB PS 16 K in the k(T) curve. The silver samples with a highest degree of perfection 
have a peak value as high as 200 W/cmK with T max «7K, but no peak at all can be found 
in the k(T) dependences of the fairly imperfect silver samples with thermal conductivity 
below 2-3 W/cmK [15]. An identical picture of the lattice-disorder influence was found 
for copper samples [28]. The thermal conductivity of the borocarbides is expected to 
be far (at least, one order of magnitude) less than that of the noble metals, since these 
latter have higher electron density and their electrons have a far greater Fermi velocity. 
The thermal conductivity of the transition metals is also far (one or even two order of 
magnitude) less than that of the noble metals. For these metals the peak in k(T) curves 
is not so sharp and is usually found only in the most perfect samples with k above 0.4-1.0 
W/cmK [15]). In some transition metals with very low thermal conductivity (Mn and V) 
the peak was not found at all [15]. It can be suggested, therefore, that the absence of a 
peak for the borocarbide k(T) curves is determined mainly by lattice defects. 

Consider now the temperature behavior of the phonon part, fttotar^e = of the ther- 
mal conductivity, presented in Figs. 3, 5, 7, and 9 by the K tota i-K e curves. According to the 
general relation (7), the k p is determined by phonon-phonon (PP), phonon-imperfection 
(PI) and phonon-electron (PE) interactions. The PP interaction should result in for- 
mation of a peak in the «p(T) dependence at T max < 0.1 Od, which, however, can be 
depressed or totaly smeared by strong enough lattice disorder. The shoulders, clearly 
seen in Figs. 3, 5, and 9 for R = Dy, Ho, and Tm in the range 35-45 K, could be indi- 
cators of the peak in the case of rather strong lattice disorder. As temperature increases 
above T max , the PP interaction can only induce a decrease in k p (T) (see Sec. Ill A). This 
is inconsistent with the temperature behavior of the n tota \-K e curves (Figs. 3, 5, 7, and 9), 
which suggests that k p rises steadily with temperature increasing up to room temperature. 
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This rise can be determined by other (PI and PE) mechanisms of phonon scattering. If 
the PI interaction is taken into account, the point defects can be assumed as the most 
important imperfections in the single-crystal borocarbides [51]. In the intermediate tem- 
perature range below 9d, the phonon wavelength decreases with increasing temperature 
approaching an interatomic distance. In this case the temperature dependence of the PI 
relaxation rate is expected to be rather weak. On the other hand, the specific heat, C p , 
of borocarbides rises significantly in this range (see discussion in Sec. IV A 1). This can 
determine (or, at least, contribute to) the rise in n p with increasing temperature. Another 
source of this increase in k v can be found in the PE interaction (see Sec. Ill C). According 
to Eq. (23), the corresponding term, n pe , grows with increasing temperature for all tem- 
peratures below Od, and this growth can be profound in the intermediate temperature 
range. It is important as well that in this range the term n pc can be comparable with the 
electronic thermal conduction (see discussion in Sec. Ill C). In this way the both, the PI 
interaction and the PE one, can contribute to the growth of the thermal conductivity of 
the borocarbides with increasing temperature. 

b) Purely superconducting borocarbides 

Up to this point the properties of magnetic superconductors (R = Dy, Ho, Er, Tm) has 
been considered. It is interesting to compare their behavior with that of the supercon- 
ducting borocarbides without magnetic order, YNi 2 B 2 C and LuNi 2 B 2 C, with Tc's 15.6 
and 16.6 K, respectively. Both of them are believed to be clean-limit type II supercon- 
ductors [32, 33, 71]. Neither Y nor Lu are magnetic and neither YNi 2 B 2 C nor LuNi 2 B 2 C 
exhibit any magnetic order [38]. No AFM correlations or local magnetic moments of 
nickel are present in YNi 2 B 2 C [72]. Inelastic neutron scattering measurements show soft 
phonon modes below T c for YNi 2 B 2 C [73] and LuNi 2 B 2 C [74]. A model phonon spectra 
was presented in Ref. [30] based on specific heat measurements [55]. A very pronounced 
low temperature softening of the low energy optical phonons and the transverse acous- 
tic modes was found. The thermal conductivity of these compounds was studied at low 
temperature in Refs. [56] (YNi 2 B 2 C) and [58, 60] (LuNi 2 B 2 C). 

The general behavior of k(T) for LuNi 2 B 2 C is shown in Fig. 11. Nearly the same 
picture is found for YNi 2 B 2 C (and not shown in this paper for this reason). The low 
temperature behavior of k(T) and p(T) is presented in Fig. 12 for both compounds. 
Figure 11 indicates that at temperature below T ~ 60 K the WF estimate for the electronic 
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FIG. 11: Temperature dependence of thermal conductivity (n tota i) for L11M2B2C. The meaning 
of K e and Ktctai-Ke is explained in the main text of the article. 

contribution, K e , is greater than the total measured thermal conductivity. 4 Obviously this 
is caused by the breakdown of the WF law in the intermediate temperature range due to 
inelastic electron-phonon scattering (see Sec. Ill B). Really, the effective Lorenz number, 
L, in this temperature range could be much less than L , used for calculation of the K e 
(see Eq.(16) and Refs. [14, 17]). The study of k(T) of a single-crystal LuNi 2 B 2 C by 
Boaknin et al. [58] also implies reduced values of L below 100 K. The deviations from 
the WF law are less for metals with increased lattice (or spin) disorder (see Sec. Ill B). 
Perhaps for this reason the breakdown of the WF law does not manifest itself so clearly 
in the magnetic borocarbides, which have higher resistivity in this temperature range 
than these of YNi 2 B 2 C and LuNi 2 B 2 C. It is interesting, however, that the K e (T) behavior 
of LuNi 2 B 2 C (Fig. 11) corresponds to one typical of electronic thermal conductivity for 
fairly good metals, and, among other things, shows a maximum at T max ~ 25 K which 
is expected at T < O.10 D (see Sec. Ill B). (0 D is about 345 K for LuNi 2 B 2 C [52]). 
This maximum (which was found for YNi 2 B 2 C as well) gives further evidence that these 
samples are well ordered. 

The superconducting transitions for LuNi 2 B 2 C and YNi 2 B 2 C are rather clearly indi- 

4 The same behavior is found for YNi 2 B 2 C for T < 80 K. 
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FIG. 12: k{T) (a) and p(T) (b) dependences for LuNi 2 B 2 C (prt/po = 24.9J and YNi 2 B 2 C 
(prt/po = 16. 5) at low temperature. The left panel (a) shows also the k(T) data for YNi 2 B 2 C 
obtained by Sera et al. [56]. 



cated in k(T) [Fig. 12(a)]. Both of the compounds show a strong enhancement in k below 
To The evident explanation is an increase in the phonon thermal conductivity due to 
reduced phonon-electron scattering as the normal electrons condense into Cooper pairs 
(see Sec. HID). The same behavior [that is a phonon-induced peak in the k(T) below 
Tc] was observed by Boaknin et al. for LuNi 2 B 2 C [58]; whereas Sera et al. [56] have 
found only a "shoulder" in k(T) below Tc for YNi 2 B 2 C which can be easily supressed by 
an external magnetic field [Sera et al. data for zero field are shown in Fig. 12(a)]. It is 
obvious that the peak magnitude and its position depends on the relative importance of 
other scattering mechanisms. In fact, the k(T) curves in Fig. 12(a) correspond to the 
case considered in Sec. Ill D when both kinds of the heat carriers (electrons and phonons) 
compete on equal terms [see Fig. 2(b)]. Taking into account Eqs. (13) and (14), the ratio 
of electronic thermal conductivities in superconducting and normal states, nl/n^, below 
Tq can be presented as 

«e/«e = n n (T)/n e , 

where n e is the electron density, n n (T) is the density of normal electrons below T c which 
decreases as temperature goes down away from To In the same way [assuming domination 
of the phonon-electron scattering and using Eq. (23)] the following relation can be written 
for the ratio of phonon thermal conductivities in the superconducting and normal states 
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below T c : 

«"/«J = n B /n n (T). 

It is seen that a decrease in n n below T c leads to continuous reduction in the ratio k%/ 
for the electronic contribution; whereas, the ratio ftp/ftp for the phonon contribution rises, 
causing an increase in the phonon thermal conductivity (see Sec. HID). The competition 
of these two mechanisms can result in the k(T) behavior shown in Fig. 12(a). 

Measurement of the thermal conductivity at very low temperatures in an applied mag- 
netic field (H cl < H < H c2 ) provides a very powerful tool for study of the symmetry of 
the superconducting gap function. At very low temperature it is somewhat easier to sort 
out the contributions to k(T) from the charge carriers and the phonons. In an ordinary 
BCS type II superconductor the normal charge carriers are localized in the vortex cores 
and can only tunnel between them in the mixed state. Thus the electronic thermal con- 
ductivity is expected to increase very slowly at low field. But as the field approaches H c2 , 
the vortices are much closer. The tunneling, then, becomes more pronounced and the 
electronic contribution increases rapidly with applied field (see general description of the 
magnetic-field effects in the thermal conductivity below Tq in Ref. [1]). If, however, there 
are nodes in the gap function, k(T,H) increases rapidly for even small fields above H c i, 
since there will be delocalized normal electrons to couple the vortex cores, even at very low 
temperature. Recent thermal conductivity measurements at low temperature for single 
crystal LuNi 2 B 2 C [60] and YNi 2 B 2 C [62] samples exhibit dramatically enhanced thermal 
conductivity at rather small applied fields. Both experiments suggest a highly anisotropic 
gap function for these borocarbides. For LuNi 2 B 2 C the data appear to suggest the pos- 
sibility of a line of nodes in the gap; whereas, the experiment with YNi 2 B 2 C has been 
interpreted in terms of point nodes along the a and b axes. These experiments certainly 
challenge the view that these materials are s-wave BCS superconductors. 

c) Purely magnetic borocarbides 

Borocarbides TbNi 2 B 2 C and GdNi 2 B 2 C are magnetic at low enough temperature 
(see Table I), but not superconducting (at least above 0.3 K). The absence of super- 
conductivity is supposed to be caused by 4/ magnetic moments [33]. The temperature 
behavior of thermal conductivity for these compounds is presented in Figs. 13 and 15. 
Expanded views of k(T), p(T) and S(T) dependences at low temperature are shown in 
Figs. 14 and 16. 

Two magnetic transitions have been found for TbNi 2 B 2 C: (1) transition (of first order) 
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FIG. 13: Temperature dependence of thermal conductivity (fttotai) f° r TbNi2B2C. The meaning 
of k c and K to tai~~K e is explained in the main text of the article. 
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FIG. 14: k(T) (a), p(T) and S(T) (b) for TbN^I^C at low temperature. The magnetic residual 
resistivity ratio, p(300K)/p(Tn) is 6.8. Using p at the lowest temperature measured, pkt/po = 
19.7 can be obtained. 
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FIG. 15: Temperature dependence of thermal conductivity (Ktotai) for GdN^E^C. The meaning 
of K e and K to tai~K e is explained in the main text of the article. 
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FIG. 16: At(T) (a), p{T) and S{T) (b) dependences for GdN^E^C at low temperature. The 
magnetic residual resistivity ratio, p(300K)/p(Tn) is 5.8. Using p at the lowest temperature 
measured, pkt/po = 29.7 can be obtained. 
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from paramagnetic to modulated AFM state below T N pa 15 K with the magnetic mo- 
ments lying along the a direction, and (2) transition (of second order) to WFM state 
below T WF =6-8 K with the magnetic moments lying along the < 100 > or < 110 > 
directions in the ab plane under applied magneic fields [38, 41, 75]. The transition to 
the modulated AFM state is clearly indicated in the k(T), and causes strong changes in 
p(T) and S(T) curves (Fig. 14). This can be explained in the same way as for DyNi 2 B 2 C 
and HoNi 2 B 2 C by an appreciable decrease in the rate of electron scattering by the spin 
disorder in response to PM-AFM transition (see discussion above). Since in this com- 
pound the electronic contribution to the thermal conductivity, K e , clearly dominates over 
the phonon one, Ktotai~^e, for the low temperature range (see Fig. 13), the measured 
thermal conductivity, K tota \, can exhibit a clear respons to this change in the electron 
scattering rate. The second-order AFM- WFM transition at T WF is not indicated in k(T) 
dependence, although weak change in p(T) and rather appreciable change in S(T) (both 
around T pa 7 K) are clearly seen (Fig. 14). The absence of some feature in k(T) at T WF 
can be explained by the fact that this transition does not strongly effect the scattering of 
charge carriers and has no influence on the phonon scattering. 

The borocarbide GdNi 2 B 2 C transforms into modulated AFM state below T N pa 20 K 
(first-order transition), with the magnetic moments lying in the b direction [76, 77]. With 
further decreasing temperature it goes into another modulated AFM state below T M pa 
13.6 K (second-order transition), with the magnetic moments tilting somewhat into the 
c direction [76, 77]. The measured K, tota i(T) dependence shows a subtle change in its 
slope at T N (Fig. 16), which is less clear than that for TbNi 2 B 2 C (Fig. 14). The weak 
feature in k(T) at T N for GdNi 2 B 2 C was found also for polycrystalline samples in Ref. 
[57]. This weak sensitivity of k(T) to first-order AFM transition can be partly explained 
by the fact that the electronic contribution to total k for this compound is much less 
than that of in TbNi 2 B 2 C (compare Figs. 13 and 15) and (especially) in DyNi 2 B 2 C and 
HoNi 2 B 2 C (see Figs. 3 and 5), where the magnetic transitions of this type produce quite 
appreciable changes in k(T) at T N (Figs. 4 and 6). The transition to another MAFM state 
below Tm pa 14 K produces only subtle changes in k(T), p{T) and S{T) of GdNi 2 B 2 C 
(Fig. 16). The explanation of this fact is the same as indicated above for TbNi 2 B 2 C, i.e., 
the transition appears to only slightly affect the scattering of charge carriers and has no 
influence on the phonon scattering. It is interesting to note that here is no indication 
of T N in the phonon contribution for R = Tb or Gd, unlike the small dip at T N in the 
phonon contribution for R = Dy shown in in Fig. 4. This suggests that the nature 
of the transitions for the former two materials is different in regards to its influence on 
the phonon contributions. An alternate explanation is that the small dip in Fig. 4 for 
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DyNi 2 B 2 C is an artifact of using the WF law in the region where p(T) is varying so 
rapidly with temperature. Any small variations in temperature between the p{T) and 
k(T) measurements would create this dip. 

d) Heavy- fermion borocarbide YbNi 2 B 2 C 

YbNi 2 B 2 C is neither superconducting (above 0.05 K [78]) nor orders magnetically 
(above 0.023 K [39]). The compound is a heavy-fermion system [39, 79, 80] with a Kondo 
temperature of T K ~ 10 K and with a Sommerfeld coefficicent of 7 ~ 530 mJ/(mole K 2 ) 
[79]. When compared to that of LuNi 2 B 2 C (7 « 19 mJ/(mole K 2 ) [52, 55]) this corre- 
sponds to an effective mass almost 30 times larger than that for the Lu-based compound 
provided that is approximately the same for all members of the series. The mag- 
netic susceptibility is anisotropic and exhibits Curie- Weiss [1/(T — #cw)] behavior at high 
temperature with a negative value of #cw, indicative that AFM correlations play a sig- 
nificant role. This behavior is consistent with the interpretation of a highly correlated 
ground state at low temperatures and crystal electric field effect at higher temperatures 
[78]. Recent neutron scattering experiments of Boothroyd et al. [81] indicate that the 
crystal electric field is significantly enhanced over that of the other family members and 
produces a temperature dependent effective Kondo interaction. These authors found a 
Kondo temperature of 25 K compared to the T K = 10 K, obtained from specific heat 
and susceptibility measurements [79], and they also suggest that at T = the compound 
might be close to a quantum critical point on the non-magnetic side. The electrical resis- 
tivity exhibits a quadratic temperature dependence below 1.5 K [79], and the ratio of the 
coefficient of the quadratic term to the 7 2 is approximately that found for UPt 3 along its 
hexagonal axis [79]. These results all suggest that hybridization between the Yb 4/ and 
conduction electron states is responsible for the suppression of superconductivity in this 
borocarbide [79] even though a simple de Gennes scaling argument [79, 80] would suggest 
that it should be superconducting at about 12 K and magnetically order at about 0.4 K. 

Single crystal alloy samples Lui_ x Yb x Ni2B 2 C have been used to explore the transition 
from superconductivity through single-impurity Kondo behavior to Kondo lattice heavy- 
fermion behavior from x = to x = 1 through transport (resistivity and thermopower) 
[82, 83] and thermodynamic (susceptibility and specific heat) measurements [83, 84]. The 
suppression of T c with Yb substitution does not scale with the de Gennes factor as would 
be expected from the Abrikosov-Gor'kov theory of magnetic pair breaking [85]. For ex- 
ample, the suppression rate of Tc at small concentrations x is 75 times higher for Yb than 
for Gd [86]. Recent studies [83, 84] indicate that the Kondo temperature varies strongly 
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FIG. 17: Temperature dependences of resistivity for YbN^E^C and LUM2B2C. The data of 
Avila et al. [88] show p(T) for a single crystal YbNi2E>2C annealed for 150 hours at 950°C 

with Yb concentration x and that the rapid suppression of T c with Yb concentration 
is consistent with Miiller-Hartmann and Zittartz theory [87] for superconducting Kondo 
systems with spin 1/2 and T K /T C ~ 10 3 . 

Although thermopower and resistivity measurements have been reported for the Lu-Yb 
borocarbide alloys, the thermal conductivity has only been reported [16] for YbNi 2 B 2 C. To 
illustrate the behavior of this system, the afr-plane resistivity, thermopower and thermal 
conductivity taken from Ref. [16] for an unannealed single crystal sample are shown in 
Figs. 17, 18, and 19, respectively. The resistivity of a single crystal sample of LuNi 2 B 2 C 
and that of the annealed YbNi 2 B 2 C sample [88] are also shown in Fig. 17 for comparison. 
At high temperature the resistivity varies approximately linearly with temperature and 
is large in magnitude compared to most metals. The magnitude of the room-temperature 
resistivity of this sample is about one-fourth of that for the polycrystalline sample in Ref. 
[80] and one-half of that for the single crystal sample in Ref. [79] . The large variations in 
the room temperature resistivity suggest unusual disorder effects. The room temperature 
resistivity for these unannealed Yb-based samples are all much larger than that for R = 
Lu, or even for the magnetic borocarbides of this family. Recent measurements [88] show 
that, although the thermodynamic properties are independent of the degree of disorder 
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FIG. 18: Temperature dependence of thermo-electric power for YbN^E^C. 

for this heavy fermion compound, the transport properties are very sensitive. Annealing 
the single crystals for 150 hours at 950°C reduced the room temperature resistivity for one 
sample from almost 100 fiQ cm to 60 pQ cm and also altered the temperature dependence 
p{T) appreciably (Fig. 17). The beginning of a Kondo minimum was observed at about 
30 K followed by the very sharp decrease in resistance due to coherence around Tk. Since 
disorder around the hybridizing Yb site can greatly change the "local" Kondo tempera- 
ture, these authors interpret the large resistivity and different temperature dependence of 
p(T) as due to a distribution of local Kondo temperatures [88]. They estimate that only a 
few percent of sites need to have Tk much greater than 10 K to explain these differences. 

The in-plane thermopower S(T) in this unannealed single crystal is negative at all 
temperatures and very large in magnitude compared to the other compounds (see Fig. 
18). It decreases linearly with decreasing temperature near room temperature and shows 
a minimum of about -44 pV/K near 70 K. At low temperature a distinct shoulder can be 
seen near 5 K in the insert. This behavior is quite different from that found in the other 
rare-earth nickel borocarbide family members [49]. These data are in general agreement 
with other measurements of S(T) for unannealed single crystals of YbNi 2 B 2 C [31, 82, 83]. 
The magnitude of thermopower, however, is also affected rather strongly by disorder, but 
the temperature dependence is not much affected. Avila et al. [88] report that annealing 
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FIG. 19: Temperature dependence of thermal conductivity (K tot3i \) for YbN^E^C. The meaning 
of k c and Ktotai-K-e is explained in the main text of the article. 

a single crystal for 150 hours at 950°C reduces the magnitude of S(T) at the minimum 
from about 40/iV/K to about 20 /iV/K and moves the minimum to somewhat higher 
temperature. The influence of disorder at lower temperatures has not been determined 
since they have measured S(T) only for T > 100 K. 

The thermal conductivity of this unannealed heavy-fermion sample is shown in Fig. 19. 
It is rather featureless except for the rapid decrease with decreasing temperature that 
starts around 100 K. Because of the large resistivity at high temperature the WF law 
is expected to be valid over much of this temperature range. The phonon contribution 
inferred from the WF law (/■Ctotai - ^) is appreciably smaller at high temperature than that 
indicated for the other borocarbide samples discussed above, but a similar increase with 
temperature at high T is observed for nearly all of the samples. As discussed earlier, 
this is most likely due to the dominance of PI and/or PE scattering. Although there is 
no k(T) data available for well annealed single crystal, one should expect singificantly 
different behavior for the electronic contribution due to the different p(T) dependence 
reported [88]. 
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B. Rutheno-Cuprates 



Coexistence of superconductivity with ferromagnetism (FM) is one of the fundamen- 
tal problems in condensed matter physics that has been studied over several decades. 
Recently, however, coexistence of these two mutually exclusive phenomena have been ob- 
served in layered cuprate systems where magnetic order first sets in at temperatures Tm as 
high as about 180 K and superconductivity as high as about 50 K. The ratio Tm/Tq ~ 4 
is much greater than that described in the intermetallic compounds discussed in sections 
I or IV A. (See Felner [89] and Lorenz, Xue and Chu [90] for a recent review.) A partic- 
ularly interesting system where coexistence of weak FM with superconductivity has been 
reported is (Ri +x Cei_ x )RuSr 2 Cu20io-5 (Ru-2122 phase), where R = Eu or Gd. The 
compounds, analogous to the Nb-2122 system [91], are related to the YBa 2 Cu 3 7 (123) 
structure with a fluorite type (Ri +x Cei_ x )02 layer separating the Cu02 double layers 
instead of the rare earth Y-layer and with Ru0 6 octahedra replacing the Cu0 4 squares 
in the CuO plane for the 123 high-Tc compound [92] . Because of this fluorite type layer, 
the successive perovskite blocks are shifted by half a diagonal of the afr-plane unit cell 
and the cell is sufficiently long to encompass two formula units per unit cell. The crystal 
structure is shown in Fig. 1 of the accompanying article by Felner [89] and Fig. 3 in the 
accompanying article by Lorenz et al. [90]. The Ru ion is pentavalent in this compound. 
Superconductivity is thought to reside primarily in the Cu02 planes, and the magnetic 
behavior appears to be associated with the Ru sites. Coexistence of superconductivity 
and weak FM was first discovered in the Ru-2122 phase for R = Eu and Gd [93, 94]. It 
was also reported for the Ru-1212 phase compound GdSr 2 RuCu 2 8 [95]. 

The only thermal conductivity measurements reported for this system are for a Ru- 
2122 sample with x = 0.5 and R = Eu which has been annealed in pure oxygen at 54 atm. 
to provide the hole doping in the Cu0 2 planes [13]. Consequently, we will focus on coexis- 
tence of weak FM and superconductivity in the compound (Eui +x Cei_ x )RuSr 2 Cu 2 Oio-5 
with x = 0.5 as evidenced by thermopower and thermal conductivity. A detailed review 
of the other properties of the Ru-2122 phase compounds is given by Felner [89] and by 
Lorenz et al. [90]. They will be only briefly summarized here before discussion of the 
thermopower and thermal conductivity. 

Hole-doping in this system of compounds can be accomplished by either adjusting 
the R 3+ /Ce 4+ ratio or the oxygen content. The parent insulator compound is assumed 
to be that for 5 = 0, x = [89]. The formal valences then associated with this com- 
pound would be Eu 3+ Ce 4+ Ru 5+ [Sr 2+ ] 2 [Cu 2+ ] 2 [0 2 "]io. Superconductivity is observed for 
0.2 < x < 0.6 with the optimally doped sample corresponding to x — 0.4. The 
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variation of Tq over this range of x is only a few Kelvin below that of the optimal 
value T c = 35 K (for x = 0.4) in the as prepared samples, much less than that ex- 
pected if all of the carriers were introduced into the Cu0 2 planes. One explanation 
that has been suggested for this is, that the holes introduced in the Cu0 2 planes by 
replacing Ce 4+ with R 3+ are partially compensated by a deficiency 5 in the oxygen 
content [89]. Neutron diffraction measurements for Gdi.3Ceo.7RuSr 2 Cu 2 Oio-<5 indicate 
8 = 0.22 with the deficiency primarily in the fluorite layer and Ru0 2 planes [96]. Re- 
cent Lni-edge XANES spectra for Gdi+^Cei^RuS^C^Oio with x = 0.9 and 0.5, es- 
tablish that Ru is pentavalent independent of the Ce concentration [89], i.e. there is 
no charge transfer to the Ru0 2 layer. In this case the formal valence count would be 
[Gd 1+:E Ce 1 _ ;c ]+ 7 - s Ru +5 [Sr 2 ] +4 [Cu 2 ] +4 - ;c + 25 [0 1 o-5]~ 20+2<5 . Thus the valence for Cu in this 
case of x = 0.3 and 5 = 0.22 would be approximately +2 as in the case of GdRuSr 2 Cu 2 8 . 
Knee et al. [96] have suggested that a hole doping mechanism arising from overlap of 
the t 2g band of Ru and the d x 2_ y 2 one of Cu, that is thought to be responsible for 
superconductivity in GdRuSr 2 Cu 2 8 , may also be responsible for superconductivity in 
Gd(Eu) 1+a .Cei_ :r RuSr 2 Cu 2 Oio-5 compounds. Hole doping can also be accomplished by 
annealing the samples at high temperature and pressure in pure oxygen. Tc has been 
raised monotonically from 34 K to as high as 49 K by annealing Eui. 5 Ce .5RuSr 2 Cu 2 Oio-5 
at 800°C in pure oxygen for 24 hours at pressures up to 150 atm. [97]. Unfortunately, 
the absolute oxygen content is difficult to determine in this compound. Consequently, 
neither Tc as a function of oxygen content nor the location of the excess oxygen in this 
compound has been systematically determined. Also, the oxygen appears to diffuse back 
out of these oxygenated samples if held in a vacuum at room temperature. 

The magnetic behavior of the Eui +:E Ce 1 _ :E RuSr 2 Cu 2 Oio-5 is also quite unusual [89]. In 
the high temperature region, Curie- Weiss paramagnetism [1/(T — 9 C w) law] is observed 
with an effective moment, g^/S(S + 1), of about 2.15 /xb for Ru and a value of 6*cw m the 
interval between 134 K and 146 K, relatively independent of the Ce concentration. This 
effective moment lies between that of the high spin state for Ru 5+ (S = 3/2, 3.9 /ib) and 
the low-spin state (S = 1/2, 1.7 //b), and the Curie- Weiss temperature #cw is positive, 
indicative of a ferromagnetic interaction. The low temperature saturation magnetization 
(gSfiB) is on iy 0.89 /^b/Ru ion compared to the value expected for the high-spin (3.0 /^b) 
and low-spin (1.0 /ib) states for pentavalent Ru. At the temperature Tm the Ru orders 
antiferromagnetically as determined by the appearance of non-linearity of in M(H). At 
a lower temperature T irr , indicated by the appearance of a remnant magnetization and 
approximate merging of the field-cooled and zero-field-cooled M(T) curves, weak FM and 
irreversibility is produced by canting of the Ru moments as a result of a Dzyaloshinsky- 



43 



Moriya type antisymmetric exchange coupling [98] between neighboring Ru moments 
that is induced by distortion of the Ru0 6 octahedrons. The sharp minimum at 77 K 
in the temperature derivative of the field-cooled DC susceptibility (xdc = M/H) shown 
in Fig. 20 coincides with this transition for a Eui. 5 Ce .5RuSr 2 Cu20io- ( 5 sample that has 
been annealed at 54 atm. in pure oxygen at 800°C for 12 hours [13]. The sharp rise in 
the derivative at 45 K coincides with the onset of superconductivity at T c . Below this 
temperature weak FM and superconductivity coexist, presumably on a microscopic scale. 
No feature indicative of T M is observed in dxDc/dT. 5 

Different possibilities for the structure of coupled inhomogeneous superconducting and 
magnetic order parameters for weak FM have been previously proposed in the literature: 
the Fulde-Ferrell-Larkin-Ovchinnikov state [99] where the superconducting order parame- 
ter develops a spatial variation with non-zero total momentum, the spiral magnetic (SM) 
state [100], the spontaneous vortex (SV) state [101], the linearly polarized (LP) state 
[102], the linearly oscillating vortex (LOV) state [103], and the spiral magnetic vortex 
(SMV) state [103]. Sonin and Felner [104] have argued that the predicted SV phase [101] 
is the most likely candidate to describe this state. 

The temperatures Tu and T irr depend rather strongly on the doping, and even on 
the method of doping. T M and T irr decrease approximately linearly with the reduction 
x in Ce content for the as prepared samples [89]. This corresponds to adding holes to 
the Cu0 2 layer, even though they may be partially compensated by an induced oxygen 
deficiency. As previously noted, there is an optimal doping for x = 0.4 which gives the 
maximum in T c . If the hole doping via Ce arises from overlap of the Ru and Cu d-bands, 
the magnetic behavior of Ru may be band magnetism instead of localized moments [96]. 
This would be consistent with the differences in the observed saturation moment and 
effective moment from those expected for Ru 5+ . It is also consistent with the absence 
of magnetic peaks in the neutron scattering measurements [89, 96]. On the other hand, 
the effects of doping by annealing in oxygen appear to be quite dependent on the Ce 
concentration. For the parent composition x — 0, there appears to be little effect on Tm 
and T irr [89]. For x = 0.5, however, both T M and T irr appear to be enhanced [97]. As 

5 Note that a totally different explanation of the magnetic behavior for this type of compound has been 
suggested by Lorenz et al. [90]. They suggest that there is a mesoscopic phase separation between 
the FM and the AFM species with superconductivity existing only in the AFM grains separated by 
nanoscale FM domains. Superconductivity would develop in the AFM domains which are Josephson 
coupled across the FM regions. The onset of the intergrain superconductivity is then expected to 
be the result of a phase-lock transition of an array of Josephson junctions rather than coexistence of 
superconductivity with ferromagnetism. 
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FIG. 20: Derivative of the field-cooled (FC) DC-susceptibility dxDc/dT [ X dc = M/H 
was recorded at H = 50 Oe] as a function of temperature T for an oxygenated 
Eui.5Ceo.5-RuSr2Cii2 0io-<5 sample (taken from Ref. [13]). The maximum value in the mag- 
nitude at about 77K, very near the point where the zero-held-cooled (ZFC) and FC curves for 
Xdc approximately appear to merge, is identified as the irreversibility temperature T| rr . The 
onset of superconductivity is identified as the sharp break in dxDc/dT at about 45K. The on- 
set of antiferromagnetic order Tm determined by the appearance of non-linearity in the M(H) 
curves at Tm = 180 K does not produce an observable feature in dxDc/dT. 

previously mentioned, Tc for this value of x was significantly enhanced by annealing in 
oxygen. Thus the effects of doping either by reducing the Ce concentration or changing 
the oxygen composition by annealing at high pressure are quite different, presumably due 
to the fact that superconductivity resides in the CuC>2 planes and the magnetic order 
resides in the Ru0 2 plane. 

The thermopower S(T) and thermal conductivity k(T) for the same sample [13] de- 
scribed in Fig. 20 is shown as a function of temperature in Fig. 21. For this sample, T M 
and T irr are 180 K and 77 K, respectively. The thermopower is positive over the entire 
temperature range, consistent with hole charge carriers, and the thermal conductivity is 
relatively small, characteristic of a pressed powder sample. There is a sharp break in 
the slope of S(T) at T = 45 K, precisely the temperature indicated in Fig. 20 for T c . 
The thermopower drops rapidly below this temperature, but it does not go to zero until 
about 29 K, as might be expected for a granular superconductor. The behavior of the 
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thermal conductivity near Tq is quite unusual, in that k(T) increases by about 30% in 
an interval of about 1 K at 48 K. Since k(T) and S(T) were measured simultaneously, 
this 3 K difference between their indicated values of Tq is probably not an experimental 
artifact. Two additional small reproducible features are seen in both S(T) and k(T), one 
at T M = 180 K, which corresponds to other indications of T M for the sample [13], and one 
at T* = 145 K, which does not coincide with any other reported magnetic or structural 
anomalies for this compound. An additional feature in k(T) is seen, a small shoulder 
with onset at T s = 13 K. Below about 5 K, the thermal conductivity begins to decrease 
rapidly. 
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FIG. 21: Thermal conductivity n{T) (a) and thermopower S(T) (b) as a function of temperature 
T for the same Eui^Ceo.^RuS^C^Oio^s sample of Fig. 20 (also taken from Ref. [13]). A 
small feature in both S(T) and k(T) at 180 K matches with the antiferromagnetic ordering 
temperature Tm, but no feature corresponding to T ur near 77K is seen in either S(T) or k(T). 
The superconducting transition temperature Tq is identified as the sharp change in the slope in 
S(T) at 45K (coincidental with that determined from dxDc/dT in Fig. 20) and as the abrupt 
jump in k{T) at 48 K. The small, reproducible feature in both k(T) and S{T) at T* = 145 K 
does not match with any yet reported magnetic or structural transitions for this compound. 

The most interesting feature is the abrupt 30% increase in k(T) at 48 K. As discussed 
in Sec. HID above, the thermal conductivity for an ordinary superconductor can either 
increase or decrease at Tq, dependent on the dominant heat carrier and scattering mecha- 
nism. If the dominant carriers are phonons and the scattering mechanism is due primarily 
to electrons, k(T) can increase below Tq as the normal electrons are frozen out. If the 
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electronic component of k dominates at Tq, an increase in k is not expected at Tq since the 
electronic part decreases below T c , as illustrated in Fig. 2. The freezing out of the normal 
electrons is generally very gradual because the conventional superconducting transition 
in the absence of a magnetic field is second order. Thus, the narrow temperature interval 
(~ 1 K) for the jump is particularly surprising. This behavior is quite distinct from any 
previously reported for high-Tc superconductors [1], which do, however, exhibit a much 
broader rise in k{T) below T c . 

Onset of magnetic order generally leads to a significant increase in k(T) due to loss of 
scattering of the electrons by spin disorder. This is particulary pronounced in the data for 
HoNi 2 B 2 C, DyNi 2 B 2 C and TmNi 2 B 2 C (though it is probably due to a new heat carrier, 
spin waves, for the latter compound) as discussed in Sec. IV A above (see Figs. 4, 6, and 
10). In contrast, there is no strong signature of magnetic order in k(T) near the AFM 
ordering temperature T N or the onset of weak FM at T WF for ErNi 2 B 2 C (see Fig. 8) which 
are both well below T c . For this rutheno-cuprate sample, however, AFM order has already 
set in at T M = 180 K and weak FM at T irr = 77 K. Based on the results described earlier 
for DyNi 2 B 2 C and TbNi 2 B 2 C (see Figs. 4 and 14), the magnetic ordering is expected to 
affect the scattering of the electrons primarily, with little effect on the phonon scattering. 
Unfortunately, resistivity was not measured for this ruthenO-cuprate sample to provide 
an estimate of the electronic contribution to k(T) through the WF law. Such an estimate 
probably would have been very inaccurate in any case, however, since it is less likely that 
this law would be valid at the higher T c values. 

The shoulder between 5 K and 13 K in k(T) is perhaps related to the similar "phonon" 
peak seen in k(T) for the Lu(Y)Ni 2 B 2 C samples discussed in Sec. IV A above (see Fig. 
12). It is quite probable that this peak for these borocarbides arises from the phonon 
channel as the normal electrons become less effective phonon scatterers at temperatures 
well below To A similar interpretation would imply a strong electron-phonon interaction 
in the rutheno-cuprate. It would also imply a very reduced density of normal electrons 
in the regions where the magnetic order parameter was strong since the shoulder was 
washed out by the introduction of vortices with their normal cores in YNi 2 B 2 C through 
application of a magnetic field [56]. The breadth of the drop to zero of S(T) at 29 K, 
well below the onset of superconductivity at 45 K, may be the result of inhomogeneity of 
the doping of the sample leading to a percolation transition or to the Josephson junction 
array model suggested by Lorenz et al. [90]. The abrupt jump in k(T) at 48 K would 
appear to be inconsistent either with a percolation transition or the Josephson junction 
array model, however. Alternatively, it may result from vortices in the superconducting 
state or another, more complicated state of the coupled superconducting and magnetic 
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order parameters. 

The thermal conductivity measurements cannot determine the detailed nature of the 
phase that leads to this abrupt jump in k near To The jump does strongly suggest 
that the transition is first order. This may be due to the influence of the internal field 
on the superconducting state, a significant change in the magnetic order to accomodate 
coexistence of both types of order, or the appearance of vortices in the superconducting 
state to accommodate weak FM. These rutheno-cuprates do appear to be quite unique 
in their behavior, however. The fact that magnetic order first sets in at a temperature 
almost four times that of T c is fascinating. A magnetic transition at 48 K to a new 
magnetic state that would permit the superconducting phase transition at 45 K would 
appear to be consistent with the abrupt jump in k and the 3 K difference in the sharp 
features near T c seen in the simultaneous S(T) and k(T) measurements. Major efforts 
are needed to grow better quality and better characterized materials so that techniques 
that probe this new state on a microscopic scale can be more readily employed. 

V. CONCLUSION 

In conclusion, we have presented here an overview of the thermal conductivity behavior 
in magnetic superconductors with particular emphasis on the rare earth nickel borocarbide 
intermetallic compounds and the 2122 rare earth rutheno-cuprates. The borocarbide 
family of intermetallics is of importance because of the availability of high quality single 
crystal samples and the large variety of phenomena - only superconductivity, coexistence of 
superconductivity and magnetic order, only magnetic order, and heavy fermion behavior - 
within the same family of compounds with identical crystal structure and similar electronic 
structure. To complete the picture of thermal conductivity in this family, a description 
of the measurements of k(T) was given for all members for which single crystal samples 
were available, not just the four magnetic superconductors. Although only limited thermal 
conductivity data is available, the oxygen rich Eui.sCeo.sRuS^Ch^Oio+s compound may 
be unique as a compound that exhibits a coupled transition to weak ferromagnetism and 
superconductivity at a surprisingly high temperature 45 K) and onset of AFM near 
200 K. 

It follows from the results presented that thermal conductivity study is a rather power- 
ful technique for revealing the fundamental properties of these complex compounds. The 
technique is especially useful for study of phase transitions. It is clear, however, that 
any single method of investigation alone cannot give a complete picture of the physical 
properties and phenomena, which are inherent in some complex compound. The ther- 
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mal conductivity studies can, however, amplify the data, obtained by other experimental 
methods (such as neutron diffraction, resistivity, specific heat, thermopower and other 
techniques) and give food for thoughts on important questions for theoretical investiga- 
tions. 
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